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ABSTRACT 

We studied the optical counterpart of the intermediate-mass black hole candidate 
HLX-1 in ESO 243-49. We used a set of Very Large Telescope imaging observations 
from 2010 November, integrated by Swift X-ray data from the same epoch. We mea- 
sured standard Vega brightnesses U = 23.89 ± 0.18 mag, B = 25.19 ± 0.30 mag, 
y = 24.79 ± 0.34 mag and i? = 24.71 ± 0.40 mag. Therefore, the source was ~ 1 mag 
fainter in each band than in a set of Hubble Space Telescope images taken a couple of 
months earlier, when the X-ray flux was a factor of 2 higher. We conclude that during 
the 2010 September observations, the optical counterpart was dominated by emission 
from an irradiated disk (which responds to the varying X-ray luminosity), rather than 
by a star cluster around the black hole (which would not change). We modelled the 
Comptonized, irradiated X-ray spectrum of the disk, and found that the optical lumi- 
nosity and colours in the 2010 November data are still consistent with emission from 
the irradiated disk, with a characteristic outer radius rout ~ 2800rin ^ IQ^^ cm and a 
reprocessing fraction w 2 x 10"'^. The optical colours are also consistent with a stellar 
population with age < 6 Myr (at solar metallicity) and mass ~ IO^Mq; this is only an 
upper limit to the mass, if there is also a significant contribution from an irradiated 
disk. We strongly rule out the presence of a young super-star-cluster, which would be 
too bright. An old globular cluster might be associated with HLX-1, as long as its 
mass < 2 X IQ^Mq for an age of 10 Gyr, but it cannot significantly contribute to the 
observed very blue and variable optical/UV emission. 
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INTRODUCTION 

•The point-like X-ray source 2XMM J011028. 1-460421 
(henceforth, HLX-1 for simplicity) is the strongest 
inter mediate-mass black hole (IMBH) cand i date known to 
date iFarrell et al.l l2009l : iDavis et all l201ll : IServillat et al] 

I2OIII ). It is seen in the sky at a distance of ~ 8" from the 
nucleus of the SO galaxy ESO 243-49 (redshift z = 0.0224, 
luminosity distance ~ 95 Mpc, distance modulus ~ 34.89 
mag; at this distance, 1" » 460 pc). Its X-ray luminos- 
ity and spectral variability (|Farrell et al.|[2009l : iGodet et al.l 

I2OO9I : IServillat et al.|[201lh are consistent with high and low 
states of an accreting BH. It shows recurrent outbursts every 
370 d (seen in 2009 August, 2010 August, 2011 August), 
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due either to some kind of disk instability, or to a periodic 
enhancement of the accretion rate. 

HLX-1 has an optical counterpart (|Soria et al.l |2010| '). 
with Hq emissio n at a red shift consistent with that of 
ESO 243-49 (IWiersema et al.ll20ia '). which strongly suggests 
a true physical association. Henceforth, we will assume that 
HLX-1 and ESO 243-49 are located at the same luminosity 
distance. The peak X-ray luminosity of HLXl is « 10 erg 
s~^, three orders of magnitude higher than typical stellar- 
mass X-ray binaries, and one order of magnitude higher than 
the most luminous ultraluminous X-ray sources (ULXs) 
found to-date. The BH mass required to be consistent with 
the Eddington limit is ~ 10* M0, and a similar value is ob- 
tained from spectral modelling of the thermal X-ray compo- 
nent, which is consisten t with emission from an accretion 
disk (jParrell et al ] I2OO9I : iDavis et all I2OIII : IServillat et al.l 

EM)- 

One of the most intriguing unsolved questions about 
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HLX-l is where it came from. There are no other lumi- 
nous X-ray sources in ESO 243-49, in the range ~ 10'^^- 
10 erg s~ ; in fact, SO galaxies do not generally contain 
extreme ULXs, which are more frequently found in spiral, 
irreg ular or interacting galaxies, with rece nt star forma- 
tion ijSwartz et al.ll2004l : PWalton et al]l201ir ). For example, 
ULXs as luminous as ~ lO'^^ erg s~^ hav e been seen in the 
Cartwheel galaxy (|Pizzolato et al.ll2010l ). representing the 
upper end of a rich population of luminous X-ray sources, 
probably high-mass X-ray binaries. Instead, HLXl appears 
to be an odd case within ESO 243-49. 

Moreover, if the BH mass estimates are correct, 
HLX-l is way too massive to have been formed from 
any stellar evolution process. Two conservative scenar- 
ios are that HLX-l may be: a) the nuclear BH of a 
disr upted dwarf satellite galaxy, accreted by ESO 243- 
49 l|King fc DehnenI bOOSl ): or, b) an IMBH formed in- 
side a massive star cluster, perhaps via stellar coa- 
lescence and merg ers (|Portegies Zwart fc McMillanI |2002| : 
iFreitag et aIll2006l V The location of HLX-l outside the disk 
plane of ESO 243-49 is consistent with both scenarios — 
in fact, massive globular clusters such as u} Cen in the 
Milky Way may themselves be the rem nants of nucleated 
dwarf galaxies (|Bekki fc FreemanI 120031 ') . IMBHs in mas- 
sive star clusters, or floati ng in the outsk i rts of gala:x- 
ies after minor mergers (|0'Learv fc Loebl |2009| ). have 
been predicted from theoretical argume nts, or specula- 
tively inferred from dynamica l modelling l|Gebhardt et al.l 



tively 1 



Liitzgendorf et al. 201l[). but their exist ence is far 



from 



proven 



van der Marel fc AndersonI |2010| ) and they 



have never been confirmed directly as luminous, accreting 
sources. Measuring the luminosity and colours of the opti- 
cal counterpart, and their long-term variability, can help us 
test those scenarios. The optical emission may come from: 
the outer regions of the BH accretion disk; a star cluster 
around the BH; a single, massive and strongly irradiated 
donor star; a combination of all three. If the optical coun- 
terpart is dominated by a star cluster, we want to constrain 
its mass and age, which may give us a clue on the IMBH 
formation process. 

In the rest of this paper, we present the results of our op- 
tical photometric study with the European Southern Obser- 
vatory (ESO)'s Very Large Telescope (VLT). We use the op- 
tical and X-ray luminosities and colours to model the spec- 
tral energy distribution, and discuss the most likely scenarios 
for the optical emission. We also compare the optical bright- 
ness in our VLT images with the brightness me asured in the 
Hubbl e Space Telescope (HST) observations of iFarrell et al.l 
l|201ll ). taken when the X-ray luminosity of the source was 
higher. 



2 VLT OBSERVATIONS 

We observed HLX-l with the Visible MultiObject Spec- 
trograph (VIMOS) mounted on the Nasmyth focus B of 
UT3 Melipal, one of the four 8-m VLT telescopes at Cerro 
Paranal. We took three images in each of the U BVRI filters. 
Our observations were carried out in service mode, and were 
split between the night of 2010 November 7 (all U, V and /, 
and two out of three B observations) and 2010 November 26 
(the last B and all R observations). All observations were 



Date 


Filter 


MJDstart 


Exp time 


2010-11-07 


U 


55508.127 


1560 s 




V 


55508.147 


375 s 




V 


55508.152 


375 s 




V 


55508.157 


375 s 




U 


55508.164 


1350 s 




U 


55508.180 


1350 s 




1 


55508.198 


165 s 




1 


55508.201 


165 s 




1 


55508.203 


165 s 




B 


55508.208 


560 s 




B 


55508.215 


560 s 


2010-11-26 


B 


55527.139 


560 s 




R 


55527.148 


285 s 




R 


55527.152 


285 s 




R 


55527.156 


285 s 



Table 1. Log of our VLT VIMOS imaging observations. 

taken in photometric conditions with an airmass < 1.2 and 
typical seeing « 0".7-0".8. Standard fields were observed 
on each night. See Table 1 for a more detailed log of our 
observations. 

The total exposure time was 4260 s in U, 1680 s in 
B, 1125 s in V, 855 s in J? and 495 s in /. In all observa- 
tions, HLX-l was located approximately at the same posi- 
tion in Quadrant 1 of the VIMOS field. We used a stan- 
dard dithering pattern for the three observations in each 
filter, to correct for pixel-to-pixel variations. For each fil- 
ter, the three images were median-combined, and standard 
calibration procedures (cosmic ray filtering, bias subtrac- 
tion, flat flelding) were applied through the ESO pipeline. 
In this paper, we focus on the Ql imaging data, although 
we also extracted and analysed the data in the other three 
VIMOS chips, which cover a significant fraction of the Ab ell 
cluster 2877. We used SExtractor (Berlin fc Arnoutll 199^ 1 to 
build a catalogue of detected objects for each of the five 
median-combined images. We t hen used SCAMP and SWarp 
iBertin et al.ll2002l : lBertinll2006t ) to read the SExtractor cata- 
logues, compute astrometric solutions, and align the images 
in the various filters. 

To convert from count rates to standard magnitudes 
(Vegamag) and fluxes at the top of the atmosphere, we used 
the standard relation 

Mag = -2.5 log (Flux [e"/s] ) + CiCol - C2 Airmass + ZP, 

where Ci is the colour coefficient suitable for a particular 
band, C2 is the extinction coefficient, and ZP is the zero- 
point (different for each filter, in the Vegamag system). We 
used a gain of 1.76e~/ADU to convert from directly mea- 
sured fluxes in ADUs/s to fluxes in e~/s. Zeropoints for our 
two nights were obtained from the pipeline-calibrated data 
files, and we also checked that they were consistent with the 
average values over the second semester of 2010, as listed on 
the ESO VIMOS quality control website Q. Average colour 
and extinction coefficients were obtained from the same web- 
site (summarized in Table 2). Not enough standard fields 
were observed in the two nights to permit us to determine 

1 www.eso.org/observing/dfo / quality /VIMOS /qc/ zeropoints.html 
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Figure 1. Top panel: combined U-band image from VLT VI- 
MOS, observed on 2010 November 7; in tiie U-band, HLX-1 is 
marginally visible by eye before background subtraction. Bottom 
panel: an example of our efforts to model ESO 243-49 with several 
GALFIT components (see text). 



Figure 2. Top row: the subsection of the VLT VIMOS U-band 
image around HLX-1. Middle row: our best-fitting thin-plate 
spline local background model. Bottom row: the background- 
model-subtracted U-band image used for determining the opti- 
cal brightness. In the left column, axis units are VIMOS pixels 
(North is up. East is left), with 1 pixel « 0".2. In the right col- 
umn, we display the same regions in 3-D format (North at the 
bottom. East to the right). 



the extinction coefficients on those specific nights; however, 
the error introduced by taking average colour and extinc- 
tion coefficients is negligible compared with other sources of 
error. Finally, we need to correct for line-of-sight Galactic 
extinction in the direction of ESO 243-49, and any additional 
extinction intrinsic to the optical counterpart of HLX-1, in 
order to determine its true luminosity and colours. We will 
discuss this in Section 4. 



3 OPTICAL DATA ANALYSIS AND RESULTS 

HLX-1 is clearly detected in emission as a point-like source 
(Figure 1), coincident with the Chandra position and with 
the location of the optical counterpart ident ified from our 
Magellan images last year (|Soriaet al.ll201(]| ). In order to 
estimate its optical flux more accurately, flrst we had to 
model and remove the contribution of the host galaxy, that 



dominates the background emission particularly at redder 
colours. To do this, we tried two alternative techniques. 

First , we tried modell i ng th e galaxy light profile with 
GALFIT (|Peng et al.ll2002l . [2OI0I ). and use the residuals to 
estimate the emission from HLX-1. The galaxy profile is 
quite complex; we tried fitting it with a combination of a 
Sersic profile plus exponential disk plus edge-on disk plus 
disky/boxy profile, and a dust mask. However, the narrow 
dust lane across the disk plane creates artifacts and rip- 
ples that badly affect the wings of the HLX-1 point spread 
function (PSF) and the determination of the background 
level in the surrounding re gion; much higher spatial reso- 
lution from HST imaging (jFarrell et al.l I2OIII ) is probably 
needed for this kind of modelling. Nonetheless, at least for 
the U and B bands, for which the galaxy halo contribu- 
tion is less important, the GALFIT residuals do show HLX- 
1 as a point-like emission source and the background sub- 



© 2011 RAS, MNRAS 000, [T]-?? 



4 R. Soria et al. 



Filter 


Zeropoint 


Colour correction 




Extinction coefficient 


Mean airmass 


U 


26.09 ±0.11 


(0.077 ±0.014) x(U - 


-B) 


0.38 ±0.04 


1.11 


B 


28.19 ±0.02 


(-0.056 ±0.005) X (B 


-V) 


0.24 ±0.01 


1.21 


V 


27.87 ± 0.02 


(-0.014 ±0.004) X (B- 


-V) 


0.12 ±0.01 


1.10 


R 


28.04 ±0.05 


(-0.102 ± 0.008) X (V 


-R) 


0.08 ±0.01 


1.20 


I 


27.64 ± 0.04 


(0.049 ± 0.006) X iV - 


'I) 


0.04 ±0.02 


1.19 



Table 2. VIMOS Ql zeropoints, colour coefficients, extinction coefficients and mean air masses used to convert from count rates to 
magnitudes. 



traction is manageable. About 3" west of HLX-1, there is 
an extended, ve ry blue source, pro bably a background star- 
forming galaxy l|FarreIl et al-lfioill ). This extended source is 
mostly responsible for the enhanc e d, asy mmetric near-UV 
emission measured by ISoria et al.l (|2010l ) in the north-east 
quadrant of ESO 243-49. After extensive tests, we concluded 
that the GALFIT method applied to our VLT data provides 
accurate results for the U and B bands, but the background- 
subtraction error is too large for the other bands. We used a 
3-pixel source extraction region, a background annulus be- 
tween 5 and 9 pixels, and applied suitable aperture correc- 
tions for the two bands (based on the PSF of isolated stars) 
to obtain the net count rates within an infinite aperture. We 
converted them to apparent magnitudes as outlined above. 
We obtained U = 23.84 ± 0.20 mag and B = 25.35 ± 0.30 
mag. 

We then used the following, more accurate technique. 
We excluded a circle of 5-pixel radius around the target 
and fitted a thin plate spline solution to the surrounding 
21 X 21 (pixel)^ area. This was done with the IDL routine 
min-curvesurf. This solution effectively follows the back- 
ground exactly outside the exclusion zone, and provides a 
smooth solution for the background inside the exclusion zone 
(Figure 2). The main caveat of this technique is that the es- 
timated background within the exclusion zone is affected by 
the values right outside the edge of the exclusion zone; in 
other words, if the exclusion circle is too small, the back- 
ground inside it is overestimated, owing to the wings of 
HLX-l's PSF extending right outside the exclusion circle. 
On the other hand, if the exclusion circle is too large, the 
fitting routine does not always provide smooth enough solu- 
tions through the whole aperture. After experimenting with 
several source and background region sizes, we decided to 
use an exclusion zone radius of 5 pixels and aperture radius 
of 3 pixels, from which we determined the net count rates 
(subsequently corrected to infinite aperture, with aperture 
corrections individually determined for each band) . We esti- 
mated that the fraction of source counts in the PSF wings, 
falling in the background region, outside the 5-pixels ex- 
clusion radius, is between 5 and 10 per cent across the 
various filters; this is negligible compared to the number 
of background counts, and does not significantly affect the 
spline fitting or contribute to the error budget. We used 
zeropoint, colour-correction coefficients and extinction co- 
efficients as before, to convert our results from count rates 
to standard (Vegamag) magnitudes, and then to AB mag- 
nitudes and flux densities (Table 3). We estimated the er- 
rors using a Monte Carlo method: we created synthetic data 
(1000 sets) based on the background model and the mea- 



sured source count rates. Modelling the simulated data in 
the same way as the real data, we obtained the photomet- 
ric errors quoted in Table 3. In summary, our flnal result 
is that U = 23.89 ± 0.18 mag, B = 25.19 ± 0.30 mag, 
V = 24.79 ± 0.34 mag, R = 24.71 ± 0.40 mag. 

HLX-1 was not detected in the I band. This was not 
surprising, given the shorter exposure time and larger unre- 
solved galactic contribution. We used a Monte Carlo tech- 
nique to estimate the upper limits for the I band flux, in 
the following manner: we took the fitted I band background 
and added a fake point-like source of increasing brightness, 
in order to see how bright a source could be "hidden" for 
a given fitted background level and measured count rate in 
the source extraction region. We found that the 2a and So- 
upper limits in the I band are / > 25.8 mag and I > 24.4 
mag, respectively (Vegamag). 

After we had completed our photometry study of the 
VLT d ata, the set of H ST images obtained and stud- 
ied by l|Farrell et al.|[201ll ) became available in the public 
archive. We used them to re-examine our source extrac- 
tion techniques, checking that there were no other point- 
like sources in the regions we chose for background model 
fitting, and testing whether HLX-1 would still be detectable 
when blurred to VLT resolution. We show (Figure 3) that 
our choice of background regions are indeed suitable. While 
HLX-1 remains easily detectable in the F390W and F555W 
even after smoothing to lower resolutions, regardless of the 
choice of region for background fitting, re-detection of the 
source in the PSF-matched F775W image (the closest match 
available to the I band) depends strongly on background lo- 
cation and modelling, which introduces other sources of sys- 
tematic error. Moreover, the source was ~ 1 mag brighter 
in every band in the //ST observations (Section 5.1). If that 
applies also to the I band, we would expect that / ~ 25.5 
mag in the VLT data; but we see from the HST F775W 
image that there are another two and possibly three faint 
sources (marked as "a", "b", "c" in the bottom panel of 
Figure 3) of comparable brightness in I (~ 25.5-26 mag), 
which may contaminate the background region. Such sources 
would make the background region brighter, and therefore 
we would oversubtract from the HLX-1 source countfl For 
these reasons, we list conservative non-detection upper lim- 
its to the I band. We stress that the other four bands with 
robust detections already provide the information needed for 
the main conclusions of our study. 



^ The same contaminating sources are also present in the other 
bands, but are a less significant problem because HLX-1 is 



brighter. 
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Filter 


(A) 


m 

(Vegamag) 


(ABmag) 


Mo 
(Vegamag) 


^obs 

(erg cm— ^ s^^ A^^) 


(erg cm-^ s-i A-i) 


(erg cm-^ s-i A— i) 


U 
B 
V 
R 
I 


3650 
4380 
5450 
6410 
7980 


23.89 ±0.18 
25.19 ±0.30 
24.79 ±0.34 
24.71 ±0.40 
> 24.4 


24.67 ±0.18 
24.99 ±0.30 
24.81 ±0.34 
24.87 ±0.40 
> 25.0 


-11.07 ±0.18 
-9.76 ±0.30 
-10.14 ±0.34 
-10.22 ± 0.40 
> -10.0 


(11.1+^;!,) X 10-19 
(5.7ii;9) X 10-19 
(4.4j:i;^) X 10-19 

< 1.7 X 10-19 


(11.8tJ;^) X 10-19 
(6.ll?;0) X 10-19 
(4.5j:};^) X 10-19 
(3.1+1;;;) ^ -^0-19 
< 1.8 X 10-19 


(14.5+^-e) X 10-19 
(r.lttr) X 10-19 
(5.2ti;9) X 10-19 
(3.4+1;^) X 10-19 
< 1.9 X 10-19 



Table 3. Summary of observed and de-reddened magnitudes and flux densities for HLX-1 during our 2010 November observations, m 
and rrtAB are the apparent magnitudes at the top of the atmosphere (Vegamag and ABmag units, respectively). Mg is the intrinsic 
absolute magnitude, corrected for a line-of-sight Galactic extinction Ay = 0.04 mag {E{B — V) = 0.013 mag) and a distance modulus 
= 34.89 mag. F"^^ is the observed top-of-the-atmosphere flux density (i.e., the flux density corresponding to m\^). is the flux density 
corrected for a line-of-sight Galactic extinction Ay = 0.04 mag (a likely lower limit to the total optical extinction). _Fi is the flux density 
corrected for an extinction Ay = 0.18 mag (a likely upper limit to the total optical extinction). The upper limits to the (non-detected) 
/ brightness and flux are at the 99% confidence level. The other error ranges are la. 



We converted from observed brightnesses in the Vega- 
mag system to top-of-the-atmosphere flux densities using 
VIMOS transmission models available onlinfl The VIMOS 
U-band filter is moderately different ( Nonin o et al.l |2009| ) 
from the Johnson U-band filter (,Bessell 199Q ) used to define 
the standard photometric system, while the other VIMOS 
filters have standard passbands. We double-checked our re- 
sults using standard Vegamag-to- ABmag conversion toolfl 
In all cases, we adopted a blackbody input spectrum at 
T = 25, 000 K for the fiux density conversion (Table 3), but 
we also calculated the flux densities in the case of a flat spec- 
trum, a T = 10, 000 K blackbody, and various stellar models 
from O to A type. The reason we adopted T — 25, 000 K is 
that the source has blue U-B and U-V colours, and does not 
appear to have a Balmer absorption jump typical of cooler 
stellar atmospheres. Anyhow, in all bands except for B, the 
fiux density at the effective wavelength of the filter (also 
known as pivot wavelength), for a given magnitude, is es- 
sentially independent of spectral shape, within ~ 1%. There 
is a small dependence on spectral shape/stellar type for the 
B band, because the Balmer jump falls at the high-frequency 
end of that passband. This introduces a systematic error of 
~ ±15% of the fiux values listed for the B band in Table 3, 
which is in any case less than the photometric uncertainty 
in the measured brightness. 

Finally, we corrected for the optical extinction, in or- 
der to determine the true emitted luminosity and colours. 
As a first approximation, we took the foreground extinc- 
tion E{B — y) = 0.013 ma.g (co rresponding to Ay ~ 0.04 
mag) from (|Schlegel et al.l 1 19981 ) . An alter native possibil- 
ity is to take the maps of Galactic HI from iKalberla et al.l 
l|2005l ). and convert A^h to a n optical extin c tion f or exam- 
ple through the relation of iGiiver fc Ozell l|2009l ). A^h = 
(2.21 ±0.09) X lO^Mv (corresponding to Ay « 0.08 mag). 
There may be additional extinction components, intrinsic 
to ESO 243-49, and in the local environment of HLX-l. For 
the ESO 243-49 component, we note that the halo of an SO 



^ http: //www. eso.org/observing/etc 

* http: / /www. stsci.edu/hst/nicmos/tools/conversion_form. html 



galaxy dominated by stellar populations older than 1 Gyr 
is not exp ected to contain much gas and dust; indeed, it 
was found (|Di Stefano et al.ll2003h that the best-fitting col- 
umn densities for the brightest X-ray sources in the Som- 
brero galaxy (an edge-on SO galaxy similar to ESO 243-49) 
are consistent or only marginally higher than the Galactic 
line-of-sight value. Moreover, optical/X-ray studies of Galac- 
tic black hole transients (fed via Roc he-lobe overfi ow, as is 
likely the case in HLX-1) have shown |Hynei[20oF) that the 
measured optical extinction is always smaller or compara- 
ble to the value that would be inferred from the values of 
A'h fitted to the X-ray spectra. This is because the local gas 
that contributes to the X-ray absorption is mostly confined 
to the X-ray emitting region, and/or is too hot to include 
a dusty component, so it does not contribute to the opti- 
cal/UV absorption. In the case of HLX-1, we found from 
the Swift/'KRT spectra that the total A'h fitted to the X- 
ray source is consistent with the foreground value, or is at 
most twic e that value (T a ble 4) . This is the same result ob- 
tained bv lServiUat et al.1 l|201lt) ( total A^h = (3 ± 1) x 10^° 
cm-^ in the thermal state), while [Codet et ali l|2009l ) adopt 
TVh « 4 X 10^° cm-^. Since we argued that the X-ray 
absorption provides an upper limit to the optical extinc- 
tion. Ay < 0.18 mag. In Table 3, we list the top-of-the- 
atmosphere fiuxes, those corrected with Ay = 0.04 mag, 
and those for Ay = 0.18 mag. The difference is small enough 
that none of our conclusions depends on the precise value of 
the extinction. 



4 X-RAY IRRADIATION 

One of our main objectives is to determine what fraction of 
the optical luminosity is due to re-emission of intercepted X- 
ray photons. Thus, we have to monitor the X-ray luminosity 
during the epoch of our VLT observations, and compare it 
with X-ray and optical luminosities at other epochs. Fortu- 
nately, HLX-1 has been the target of more than 100 Swift 
X-ray Telescope (XRT) observations since 2008 October; see 
NASA's HEASARC data archive for a detailed logbo ok. We 
used the on-line XRT data product generator () Evans et al.l 



© 2011 RAS, MNRAS 000, [T]-?? 



6 R. Soria et al. 



HSTAVFC3 F390W 



4".3 





Figure 3. Top panel; _ffST/WFC3 image with the F390W fil- 
ter, taken on 2010 September 23, with the outline of the region 
chosen for background modelling in our VLT study. North is up. 
East to the left. Middle panel: same as in the top panel, for the 
F555W filter image. Bottom panel: same, for the F775W filter 
image. In all panels, the box size is 4". 3 X 4". 3 and the ra- 
dius of the exclusion circle is 1". Cosmic rays inside or near 
the extraction regions are marked with yellow crosses. The three 
faint sources marked with "a" , "b" , "c" in the F775W image may 
affect background subtraction in the VLT I band. 



l2007l . l2009h to extract light curves and spectra (including 
background and ancillary response files); we selected grade 
0-12 events. We downloaded the suitable spectral response 
file for single and double events in photon-counting mode 
from the latest Swift Calibration Database. 

In 2010 November, HLX-1 was in the decline phase af- 
ter the 2010 X-ray outburst. The X-ray luminosity did not 
change significantly between 2010 November 7 and 2010 
November 26 (Figure 4), and the timescale for the expo- 
nential decay is much longer than 19 days. Thus, we can 
assume that the irradiation correction to the R band (ob- 
served on November 26) should not be significantly different 
from that affecting the optical bands on November 7. The in- 
dividual Swift/X'RT observations around the two epochs are 
too short to allow for meaningful spectral analysis. There- 
fore, to characterize the X-ray spectrum, we coadded all the 
observations in the middle part of the outburst decline, be- 
tween 2010 September 30 and 2010 December 11, a period in 
which the X-ray source was almost on a plateau. Moreover, 
we coadded all Swift/XRr observations from 2008 October 
to 2011 August in which the 0.3-10 keV count rate was be- 
tween 0.01 and 0.02 ct s~^. We then fitted bot h co-added 
spectra with XSPEC version 12 (|Arnaud|[l996l). using the 
Com ptonized, irradiated disk model diskir iGierlmski et al.l 
l2008l ). 

An advantage of the diskir model over widely used two- 
component models such as diskhb plus power-law is that the 
Comptonized component is physically truncated at low en- 
ergies, and the disk spectrum can easily be extrapolated to 
the optical bands, and compared with the flux from the op- 
tical counterpart. We fixed the fraction /in of luminosity in 
the Compton tail that is thermalized in the inner disk to 
0.1 (|Gierliriski et al.ll2008t ). The other irradiation parameter 
in the model, /out, is the fraction of bolometric flux that is 
thermalized in the outer disk. The value of /out and of the 
outer disk radius rout are essentially unconstrained when 
we fit the X-ray data, especially for a cool disk with peak 
temperature < 0.2 keV, but will become important for the 
optical extrapolation. 

We find that the long- and short-baseline Suii/i/XRT 
spectra are essentially identical in shape and normalization 
(Figure 5). The exposure time for the long spectrum was 
1.2 X 10^ s; the exposure time for the short spectrum (which 
is simply a subset of the longer-baseline data) was 4.7 x 10* 
s. The best-fitting parameters for the long-baseline spec- 
trum are listed in Table 4. This result is consistent with the 
expected spectral behaviour o f HLX-1 at intermediate lumi- 
nosities (|Servillat et al.ll2011^ . with a soft thermal compo- 
nent that is still dominant, and an emerging power-law that 
becomes harder during the decline. The emitted luminosity 
Lo.3_iokcV ~ 5.5 X 10*^ erg , a factor of 2 lower than at 
the outburst peak. The disk normalization implies an inner 
radius rin « 6 x 10'* km (for an approximately face-on disk), 
indicative of a BH mass ~ lO^Mp and consistent with pre - 
vious X-ray studies (|Davis et al.ll201ll : IServiUat et al.ll201ll '). 



5 OPTICAL/X-RAY COMPARISON 

After measuring the optical brightness in the VLT dataset, 
and modelling the temperature and luminosity of the X-ray- 
emitting accretion disk, we shall now discuss whether the 
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Figure 4. Sm/t/XRT lightcurve in tiie 0.3-10 keV band, be- 
tween 2008 October 24 and 2011 November 28. The two epoclis 
of our VLT observations are indicated by vertical green lines; the 
epoch of the HST observations is also plotted. We integrated to 
Swift/XRT spectrum from the epochs plotted in red, to determine 
the irradiation conditions around the time of our VLT observa- 
tions. 
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Figure 5. Swift/XRT spectrum modelled with the Comptonized, 
irradiated disk model diskir in XSPEC. The blue spectrum in- 
cludes all the epochs when the XRT count rate was between 0.01 
and 0.02 ct . The red datapoints are a subset of the blue data- 
points, for the epochs plotted in red in Figure 4. The residuals 
in the bottom panel refer to the long-baseline spectrum, which 
has better signal-to- noise and provides stronger constraints. 



optical emission is more likely to come from the outer part of 
the irradiated disk, from a stellar population or from a single 
donor star. A com parison with the op tical brightness in the 
//ST observations fFa rrell et al ]|201ll ') provides a crucial test 
to distinguish between those models. 



Table 4. Best-fitting spectral parameters for the coadded 
Swift/XRT observations during intermediate-luminosity states. 
Units: TVjj Gal a^nd AfH,int ^.re in units of lO'^" cm~'^; feTin and 
kTe in keV; ri^j in units of the inner-disk radius rj^; /o.3— lo in 
10"!'^ erg cm~2 s"^; Lo.3-iO in 10*^ erg s^^. The disk normal- 
ization Af^iak = [(''in/km) / (D/lOkpc)] ^ COS S, where cos9 is the 
viewing angle of the disk. Errors indicate the 90% confidence in- 
terval for each parameter of interest. 

5.1 Optical decline between the HST and VLT 
observations 

The X-ray luminosity, as determined from the Swift/XKT 
count rate and spectra, was a factor of « 2 higher during 
the 2010 September 23 HST observations than during our 
VLT observations, taken at a later stage of the 2010 out- 
burst decline. We showed (Table 3) that the optical/near-UV 
brightness of HLX-1 was « 24.7-25.0 mag (ABmag system) 
in out VL T data; in the HS T observations, it was « 23.8- 
24.0 mag (|Farrell et all 1201 ll fl. Ahhough the VIMOS filter 
bands are slightly different from the HST bands, this dis- 
crepancy in flux density occurs in every band and cannot 
be explained with colour transformations; as we are refer- 
ring to the top-of-the-atmosphere flux density, this is also 
independent of the different choices of line-of-sight and in- 
trinsic extin ction correct ions. Interpolating the flux densi- 
ties given bv lFarrell et al. (2011) across our observed bands, 
we determine that the VLT flux was fainter than during 
the HST observations, with 3cr signiflcance in each of the 
U,B,V and I bands, and 2.50- significance in R. Thus, we 
conclude that the source had significantly dimmed, by « 1 
mag, in the near-UV/optical/near-IR bands. To test possi- 

^ We also independently checked that HLX-1 was significantly 
brighter in the HST images, by doing relative photometry with a 
dozen other sources detected in both fields. 
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ble further evidence of variability in the optical c ounterpart, 
we al so re-examined our Magellan observations (|Soria et alj 
taken a few days after the 2009 August outburst 
peak, when the X-ray flux was, again, a factor of 2 higher 
than in 2010 November. The Magellan R-band brightness 
mAB.ij ~ 24.0 ± 0.3 mag is ~ 0.9 mag higher than in our 
2010 VLT observations, with a 2a significance. (The larger 
error in the Magellan V-band measurement, which we now 
re-estimate as mAB,v « 24.5 ± 0.5 mag prevents firm con- 
clusions on that band.) 

On its own, our discovery of long-term optical variabil- 
ity correlated with the X-ray luminosity already rules out 
the possibility that the optical/UV counterpart was domi- 
nated by emission from an old globular cluster or a young 
su per-star-cluster dur ing the i/ST observations, as discussed 
in lFarrell et al.l (|201ll ). Hence, we suggest that the dominant 
optical/UV contribution near the outburst peak was from an 
irradiated accretion disk. Any additional contribution from 
a (constant-luminosity) stellar population must be fainter 
than the luminosity determined from our VLT observations. 

In fact, even in our VLT data, the accretion disk is likely 
to give a significant contribution to the optical/UV fiux, as 
the X-ray luminosity was still 20 times higher than in the 
low state, and the soft X-ray flux was still dominated by 
(partly Comptonized) disk emission (Table 4). In the rest 
of this Section, we are going to discuss and constrain the 
possible contribution of disk and stellar populations in our 
VLT data. 



5.2 Irradiated disk 



extending to rout ~ 2800rin ~ (1.7 x 10^^)/ cos ^ cm, and 
with an intercepted and re-radiated fraction of X-ray pho- 
tons /out ~ 2 x 10~^, a very plausible value for irradiated 
disks in binary systems. Reprocessing fractions of ~ a few 
times 10~'^ are suggested both by theoretical modelling (e.g. . 



The simplest scenario is a standard (jShakura fc SunvaevI 
1 19731 ) non-irradiated disk, with inner temperature and nor- 
malization derived from X-ray fitting (Table 4), extrapolated 
to an infinite radius (for practical purposes, rout = lO^rin). 
This is plotted as a solid black line in all panels of Figure 6; 
in those plots. X-ray absorption has been removed from the 
disk model, and the optical datapoints have been corrected 
for a line-of-sight reddening E{B -V) = 0.013 mafl In the 
optical band, the predicted non- irradiated disk emission ap- 
proximately reproduces the spectral slope between U and R 
but falls a factor of 3 below the observed brightnesses. 

Clearly, irradiation of the outer disk is needed to in- 
crease the optical fiux of the model without changing the 
X-ray luminosity. First, we assumed an irradiation fraction 
/out = 10"'^ while keeping a large outer radius (rout = 
lO^rin). The predicted optical fiux (solid green line in Fig- 
ure 6) is now roughly consistent with the VLT data from U 
to R (although the U-V colour is too red), but is much too 
bright for the I band. As a general rule, the optical emis- 
sion of a large irradiated disk is redder than the emission 
of a non-irradiated disk. However, our source is significantly 
bluer. 

Truncating (or shading) the outer edge of the disk re- 
moves some of the red emission from the outer regions, and 
makes the optical slope steeper. We find that the optical 
data are well fitted (solid blue line in Figure 6) with emis- 
sion from the same accretion disk seen in the X-ray band. 



Vrtilek et al lll99oi : 1de Jong et al.|[l996l : iKing fc Ritteilll998l : 
Dubus et al.ll 19991 ) and by obs ervations of a ccretion disks i n 



Galactic BHs fe.g.. lHvnes"eral...2002 : Gierliiiski et al.ll2009l '). 
In summary, considering that an irradiated disk was the 
dominant optical/UV emitter in 2010 September (Section 
5.1), and that the X-ray flux from the inner disk was still 
relatively high in 2010 November, we take the irradiated disk 
scenario as the most likely interpretation also for the optical 
emission in our VLT data. 



5.3 Stellar population 

An alternative possibility is that most of the optical emis- 
sion in the VLT data comes from a massive star cluster sur- 
rounding the BH. For example, this is what we expect if the 
IMBH was the nucleus of a recently accreted satellite dwarf, 
or was formed from core collapse in a massive star clus- 
ter. We ran inst antaneous star formation simulations with 
STAR BURST99 (|Leitherer et al.lll999l : IVazguez fc Leithereij 
to determine the expected brightness and colours of 
star clusters of various masses and ages (we assumed solar 
metallicity for simplicity, and a Kroupa initial mass func- 
tion). Here, we focus on a few representative cases. 

We find that a young super-star-cluster (for example, 
the solid cyan line in the top panel of Figure 6 represents a 
cluster with M = l(f Mq and age of 10 Myr) is inconsistent 
with the VLT optical data: it would have blue colours and 
small or no Balmer jump, consistent with the observed data, 
but would be muc h too bright. This ru les out one of the sce- 
narios presented in lFarrell et al.l l|201ll ). A young star cluster 
with much lower mass can match the observed optical lumi- 
nosity and colours; the solid cyan line in the bottom panel 
of Figure 6 represents a cluster with M — 1.4 x 10** and 
age of 5 Myr. In fact, this is only an upper limit to the young 
cluster mass, if there is still a significant disk contribution 
(Section 5.2). A stellar population age < 6 Myr is required 
in order to match the optical colours, at solar metallicity — 
modulo the age-metallicity degeneracjQ. However, the phys- 
ical interpretation of this scenario is problematic. A IG^Mq 
cluster is too small to harbour an IMBH formed from 
runaway core colla pse (|Portegies Zwart fc McMillanI [20o3 : 
iFreitag et al.ll2006|). especially one as massive as required 
for HLX-1 ( Davis et al.l 12011 ). similar to the mass of the 
cluster itself. Its very young age makes it unlikely that the 
cluster could be much more massive at birth. Moreover, the 
HST near-UV images show no other signs of ongoing or very 
recent, clustered star formation in the region surrounding 
HLX-1, which appears like a one-off object. For all these 
reasons, we conclude that the young-star-cluster scenario is 
much less likely than the irradiated disk scenario. 

Old globular clusters are much fainter for a given stellar 
mass. A cluster with M — IQ^Mq and age of 1 Gyr would be 



^ Using our upper limit E{B — V) = 0.06 mag docs not change 
any of tlie arguments in this Section. 



^ In simple terms, the age-metallicity degeneracy means that 
colours and luminosity are the same for a stellar population wit h 
half the metal abundance and three times the age iworthevll994) . 
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Figure 6. Top panel: optical/UV flux densities from our VLT 
data (plotted in red), compared with the model flux density pre- 
dicted from various types of accretion disks and stellar popula- 
tions. At high energies, the disk parameters are fixed from the fit 
to the Swift/XRT data. Black line: standard, non-irradiated disk 
spectrum; green line: irradiated disk; blue line: irradiated, trun- 
cated disk; cyan line: spectrum of a young, massive star cluster 
magenta line: spectrum of an intermediate-age stellar population: 
orange line: spectrum of an old globular cluster. Middle panel 
zoomed-in view in the optical/UV band, plotted in linear scale 
the three disk spectra and the two old/intermediate-age stellar 
populations are as in the top panel. Bottom panel: the spectrum 
of a low-mass, very young stellar population (cyan line) fits the 
data as well as the irradiated, truncated disk model (blue line) 
but may not have a simple physical interpretation; an very old, 
massive globular cluster may also be associated with HLX-1, pro- 
vided that its mass is < 2 X 10^ Mq (magenta line). 
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Figure 7. Optical/UV flux densities from our VLT data (plotted 
in red), compared with a model consisting of a standard non- 
irradiated disk (dotted blue line) plus a single blackbody with 



T = 20, 000 K and characteristic size 
(dotted green line). 
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marginally consistent with the observed B, V, R brightnesses 
(solid orange line in Figure 6, middle panel), but totally in- 
consistent with the observed colour U — B ~ —1.3 mag, 
which is typical of OB stars, and with the upper limit to 
the I band. Intermediate-age (« 50-100 Myr) clusters with 
masses ~ W^Mq are a slightly better match to the B, V, 
R optical data but remain inconsistent with U and I (solid 
magenta line in Figure 6, middle panel). The strongest con- 
straint against intermediate-age and old stellar populations 
comes from the apparent lack of a Balmer jump between B 
and U, in our VLT observations; this implies that Hydrogen 
is mostly ionized at the photosphere of the optical emitter. 
If the emitter is a stellar population, the lack of a Balmer 
jump alone strongly suggests an age < 50 Myr. We conclude 
that the optical counterpart of HLX-1 cannot be dominated 
by a massive old or intermediate-age globular cluster. 

More complex and/or contrived theoretical scenarios 
can also be proposed. For example, the IMBH could be lo- 
cated in an old, faint globular cluster (perhaps the parent 
cluster that formed the IMBH via core collapse in earlier 
epochs) with an age > 1 Gyr and current mass < a few 
xW^Mq (thus, satisfying the VLT luminosity constraints), 
which has recently undergone a new burst of star forma- 
tion, that produced a ~ IO'^Mq young stellar population 
(including perhaps the donor star to the IMBH). Or, the 
blue optical emission may come entirely from the IMBH 
disk, inside an old globular cluster, and the mass donor may 
be an evolved low-mass star. For a globular cluster with an 
age of 10 Gyr, the upper mass limit consistent with the VLT 
data is « 2 X 10^ (Figure 6, bottom panel). Such an old 
cluster would dominate the near-IR emission but would give 
a negligible contribution to the optical/UV bands, compared 
with the irradiated disk and/or young stellar population. We 
do not have enough observational data at the moment to 
test such multiple-component scenarios. New observations 
of HLX-1 during a much lower X-ray state, and in the IR, 
will provide conclusive tests between the irradiated disk and 
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star cluster models, and may provide a more quantitative es- 
timate or upper limit to the optical contribution of a stellar 
population around the BH. 

Finally, we may want to consider the possibility that 
the optical emission comes from the surface of an individual 
donor star. The inferred absolute magnitude Mv ~ Mb ~ 
— 10 mag is very high but not beyond the limit of exist- 
ing stars, for example blue hypergiants and Luminous Blue 
Variables (LBVs). Stars in this luminosity range include 
Cyg OB2-304, the most lumino us star known to -date in the 
Galaxy {My ~ —10.6 mag: Massev et al.ll200j ): r; Carinae 
was also at My ^ — 10 ma g in the ear ly 19th century, before 
the 1843 Great Eruption (jFrewll2004l '). 

As an extreme case, we modeled a standard non- 
irradiated disk (based on the X-ray parameters, as in Sec- 
tion 5.1) with an additional single-temperature blackbody 
component. We find that for Tbb ~ 20, 000 K and a stellar 
radius r, ~ 8.9 x 10^^ cm ~ ISORq, the resulting spectrum 
also matches the observed optical data (Figure 7) as well as 
the irradiated disk spectrum. We cannot distinguish between 
the two cases with the scant optical data at hand. What this 
means is that most of the optical emission may come from a 
region with characteristic size ~ 10^^ cm and characteristic 
temperature ~ 20, 000 K: either the outer radius of the ir- 
radiated (part of the) disk, or the irradiated face of a large, 
massive donor star. A scenario where the optical luminos- 
ity is dominated by a single LBV or blue hypergiant donor 
star implies the presence of ellipsoidal variations, due to the 
tidally-distorted shape of the star, and orbital modulations, 
due to the variable fraction of the irradiated face of the star 
visible to us during the orbit. The single-donor scenario may 
perhaps explai n the outburst peri od ~ 370 d seen in the X- 
ray lightcurve (|Lasota et al.|[201ll ). if the orbit is elongated 
and Roche-lobe mass transfer occurs mostly at periastron, 
with the formation of a disk. However, it is very unlikely that 
the optical/UV counterpart is a single donor star, because 
there are no other similar objects or other evidence of ex- 
ceptional star formation in the surrounding regions, within 
at least 1 kpc. Another issue to consider is that such mas- 
sive stars are shrouded in dense winds, in apparent contrast 
with the low column density measured from the X-ray spec- 
tra; on the other hand, the wind may be completely ionized 
because of the X-ray emission from the BH. Deeper opti- 
cal spectra of the optical counterpart may help test between 
those scenarios. 



6 CONCLUSIONS 

We observed the optical counterpart of HLX-1 with VLT 
VIMOS, in 2010 November, when the X-ray source was in an 
intermediate state, declining after the 2010 August outburst. 
We found the following brightness values in the standard 
Vegamag system: U = 23.89 ± 0.18 mag; B = 25.19 ± 0.30 
mag; V = 24.79 ± 0.34 mag; R = 24.71 ± 0.40 mag. In the 
ABmag system, this corresponds to rriAB.u = 24.67 ± 0.18 
mag, mAB,B = 24.99 ±0.30 mag, mAB,v = 24.81 ±0.34 mag, 
wiAB,R = 24.87 ± 0.40 mag. We do not detect the source in 
the I band (mAB,i > 25.0 mag); this was expected because 
the I band had the shortest exposure time and the largest 
galactic contamination. 

We found that the optical/UV brightnesses in every 



band are ~ 1 mag fainter than in the HST observations from 
2010 September, at a time when the X-r ay fiux was a factor 
of 2 higher. This rules out the scenario (|Farrell et al.|[201lh 
that the optical source seen in the HST observations was 
dominated by a massive star cluster, because in that case 
its brightness would not have changed. It strongly suggests 
that the optical counterpart was dominated by an irradiated 
accretion disk. 

From the VLT data alone, we cannot tell whether the 
optical emission was still dominated by an irradiated disk, 
or by a stellar population: optical observations at even lower 
X-ray fiux levels will test these alternatives. We can, how- 
ever, constrain the maximum stellar mass and age of the 
stellar population that can be consistent with the observed 
VLT luminosities and colours. We found that they are con- 
sistent with a young stellar population of mass ^ 
and age < 10 Myr. This immediately rules out a young 
super-star-cluster, which would have similar blue colours 
and lack of a Balmer jump, but would be much too lumi- 
nous. If there is an old globular cluster around HLX-1, as 
previously speculated ((Soria et al.ii201Q) . it may dominate 
the near-IR bands, but it cannot give a significant contribu- 
tion to the observed optical/UV flux discussed in this work. 
We estimate that a 10-Gyr-old globular cluster must have a 
mass < 2 X IO^Mq to be consistent with the VLT data. By 
comparison with the current models and observational mass 
limi ts on IMBH candidates in old gl obular clusters (Table 
1 in Ivan der Marel fc Andersonll2010l ). we argue that a host 
cluster mass « 2 x 10® Mq may be just high enough at least 
to allow for the possibility of a central IMB H with a mass 
~ lO '^ Mp), as suggested by X-ray modelling (|Servillat et al.l 
l201ll : lDavis et al.ll201ll ). 

In conclusion, we attribute most of the optical/UV 
emission to an irradiated disk, as the simplest scenario con- 
sistent with all the data at hand. The more complex scenario 
of a very old globular cluster with a small burst of recent star 
formation is also acceptable, but requires more observational 
data to be tested. We modelled the X-ray spectrum around 
the time of the VLT observations (when Lo.s-io ~ 5.5 x lO*'^ 
erg s"'^), and during epochs of similar levels of X-ray lu- 
minosity in previous years. We showed that an irradiated 
disk can explain the observed optical luminosity and colours, 
with a plausible reprocessing fraction « 2 x 10"'^, typical of 
irradiated disks in X-ray binaries. Regardless of any details 
on disk structure and irradiation, if the optical emission is 
thermal, during the VLT observations it comes from a region 
with characteristic size ~ 10^'' cm and characteristic tem- 
perature ~ 20, 000 K; this could be the size and temperature 
of the irradiated outer disk. 
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